import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
from math import sin, cos, pi

x = np.linspace(1, 2, 100)
y = np.linspace(1, 2, 100)
X, Y = np.meshgrid(x, y)
Z = np.sin(X)*np.sin(Y)

fig = plt.figure()
plt.xlabel('x')
plt.ylabel('y')
ax = Axes3D(fig)
ax.plot_surface(X, Y, Z)
plt.show()
